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Abstract. We previously obtained constraints on the viewing geometries of 6 Fermi LAT 
pulsars using a multiwavelength approach ,1 . To obtain these constraints we compared the 
observed radio and 7 -ray light curves (LCs) for those 6 pulsars by eye to LCs predicted 
by geometric models detailing the location and extent of emission regions in a pulsar 
magnetosphere. As a precursor to obtaining these constraints, a parameter study was conducted 
to reinforce our qualitative understanding of how the underlying model parameters effect the 
LCs produced by the geometric models. Extracting useful trends from the 7 -ray model LCs 
proved difficult though due to the increased complexity of the geometric models for the 7 -ray 
emission relative to those for the radio emission. In this paper we explore a second approach to 
investigating the interplay between the model parameters and the LC atlas. This approach does 
not attempt to understand how the set of model parameters influences the LC shapes directly, 
but rather, more fundamentally, investigates how the set of model parameters effects the sky 
maps from which the latter are extracted. This allows us to also recognise structure within the 
atlas itself, as we are now able to attribute certain features of the LCs to specific features on 
the sky map, meaning that we not only understand how the structure of single LCs come about, 
but also how their structure changes as we move through the geometric solution space. 


1. Introduction 

The Fermi Large Area Telescope (LAT) is the primary instrument on the Fermi Gamma-ray 
Space Telescope mission which was launched on 11 June 2008. Fermi LAT is an imaging, wide 
field-of-view, high-energy 7 -ray telescope covering the energy range from 20 MeV to 300 GeV. 
The Fermi LAT continuously scans the sky in an ever improving all-sky survey, taking advantage 
of it’s large field-of-view and high sensitivity to scan the entire sky every three hours. One of 
the important products to date of this all-sky survey is the first Fermi LAT catalogue of 7 -ray 
pulsars \2\ produced using the first six months of Fermi LAT data, and increasing the number 
of known 7 -ray pulsars from at least 6 [3] to 46. The second Fermi LAT 7 -ray pulsar catalogue 
is due to be published soon, and will push the number of known 7 -ray pulsars to over 100. This 
dramatic increase has opened up a whole new population of pulsars, with exciting possibilities 
for multiwavelength studies [ 3 ] . 

We have previously obtained fits to 6 Fermi LAT pulsar LCs using both radio and 7 -ray data 
to concurrently fit radio and 7 -ray geometric pulsar model LC predictions [5j. 


1.1. The geometric pulsar models 

We use an idealized picture of the pulsar system, wherein the magnetic field has a retarded 
dipole structure [6j and the 7 -ray emission originates in regions of the magnetosphere (referred 
to as ‘gaps’) where the local charge density is sufficiently lower than the Goldreich-Julian charge 
density |7]. These gaps facilitate particle acceleration and radiation. We assume that there are 
constant-emissivity emission layers embedded within the gaps in the pulsar’s corotating frame. 
The location and geometry of these emission layers determine the shape of the 7 -ray LCs, and 
there exist multiple models describing their geometry. 

Two such models for the geometry of the 7 -ray emission regions are the outer gap (OG, [HI) 
and two-pole caustic (TPC, [9]) models. In both the OG and TPC models emissions are produced 
by accelerated charged particles moving within narrow gaps along the last open magnetic field 
lines (the field lines which close at the light cylinder). In the OG model, radiation originates 
above the null charge surface (where the Goldreich-Julian charge density changes sign) and 
interior to (closer to the magnetic axis) the last open magnetic field lines. The TPC gap 
starts at the stellar surface and extends along the last open field lines up to near the light 
cylinder, where the corotation speed approaches the speed of light. The special relativistic 
effects of aberration and time-of-flight delays, which become important in regions far from the 
stellar surface (especially near the light cylinder), together with the curvature of the magnetic 
field lines, cause the radiation to form caustics (accumulated emission in narrow phase bands). 
These caustics are detected as peaks in the observed 7 -ray LCs B HD|. 

1.2. Phaseplots and LCs 

The tool we use to visualise the emission of the pulsars is called a phaseplot, and is an 
equirectangular projection skymap of the emission of the pulsar with the phase at which the 
emission is observed, </>, on the horizontal axis, and £ on the vertical axis. For each point in the 
(</>, () space it gives the relative intensity of the emission per solid angle which would be observed 
if the line-of-sight to the observer would pass through that point. On this phaseplot the pulsar 
makes one full rotation around its rotation axis, and the equator of the NS is at £ = 90°. To find 
out what an observer would see during the course of the period we can obtain an LC from this 
phaseplot by taking a cut at a constant ( through the phaseplot. Fig. [2] shows an example of a 
phaseplot (for the TPC case) for a pulsar with an inclination angle of a = 50°, accompanied by 
an LC extracted from it using such a cut (designated by the yellow line through the phaseplot). 
Both the phaseplot and LC are independently normalised so their respective maxima are 1. 

2. Obtaining LC fits 

To obtain fits to the observed pulsar LCs we need to consider the entire domain of a and £ at the 
very least, keeping the other parameters fixed. To most effectively do this we generate ‘atlases’, 
consisting of a set of model LCs covering the (a,£) domain at a chosen resolution, i.e., a grid 
of LCs sampling the two-dimensional parameter space evenly. The data are then compared to 
this atlas (for our purposes this is done by eye) and the best-fit (cq Q is then inferred. 

This fitting process has proven to be effective in a multi wavelength arena with concurrent 
fitting of both radio and 7 -ray light curves yielding well-constrained best-fit parameters for 
a number of radio-loud 7 -ray pulsars jT]. In the case of by-eye fitting the effectiveness and 
usefulness of this process is quite dependent on the level of the qualitative knowledge available 
about the behaviour of the produced atlas of LCs when model parameters other than a and £ 
are varied. This knowledge can be best garnered by performing parameter studies on the models 
focusing on the relevant parameters. 


2.1. Performing parameter studies: radio models vs. 7 -ray models 

To perform parameter studies using these phaseplots and LCs it is convenient to investigate the 
interplay between the parameters and the LC atlas. In the case of the simpler conal radio models 
this approach proves very effective and a lot can be learned from doing such parameter studies. 
The 7 -ray models yield considerably more complex phaseplots, and as such, more convoluted 
LC atlases. This means that to obtain useful knowledge from these parameter studies we need 
to adjust our approach accordingly. 



Figure 1 . The process whereby fits to observed LCs can be obtained from implemented 
geometric models. The red arrow denotes the methedology of understanding the structure 
of the atlas in light of only the structure of the LCs. The blue lines (solid and dashed combined) 
denote how we propose this understanding of LC structure could be attained more effectively. 

Fig- [U shows a schematic representation of the process used to obtain fits to the radio and 
7 -ray LCs. The red arrow illustrates the established methodology of understanding the influence 
the parameters have on the LC atlas, by examining the LCs directly. The blue arrows illustrate 
the proposed methodology wherein the influence of the model parameters on the LC atlas is 
understood by examining the effects of the model parameters on the structure of the phaseplots 
from which the LC atlas is generated. 

3. Re-examining the LCs and phaseplots 

Fig-i shows how the different features on the phaseplot can be associated with the peaks seen 
on the LC. In this case (the TPC model with a = 50° and £ = 70°) the two main peaks 
on the LC are produced by cutting the caustic (marked C on the phaseplot) above its lowest 
point (around £ = 50°), thus intersecting it twice. The earlier (at lower <f>) of these two peaks 
also includes emission from the so-called ‘overlap region’ (where the directions of many of the 
rotationally distorted magnetic field lines are such that their emission overlaps; marked O on 
the phaseplot). The LC furthermore includes four smaller peaks, one preceding the earliest large 
peak, one between the two large peaks, and two following the later large peak. The earliest of 
the four is associated with the so-called ‘dynamical caustic’ (marked D on the phaseplot), and 
is the dominant peak for LCs produced cutting below (at smaller Q the lowest point of the 
caustic. The second and fourth peaks of the four smaller peaks are associated with the caustic 
itself, respectively being produced by very high altitude emission (where the magnetic field is 
most distorted), and by the small tail of the caustic (both encircled). The third of the four 
smaller peaks is produced by the line on the phaseplot associated with the so-called ‘notch’ in 
the rim of the neutron star’s polar cap at the surface. 

Thus, we see that the parts of an extracted LC can be understood in terms of the features 
on the phaseplot which generate them. This means that the hard work of understanding how 
the LCs change under changes in the underlying parameter set, shifts to having to understand 
how the structure of the phaseplots change under changes in the underlying parameters. There 
are mainly 5 features on the phaseplot we need to account for, namely: 
































Figure 2. A phaseplot generated by the TPC model with a = 50° and cut at ( = 70°. C 
denotes the main caustic, D denotes the secondary caustic, P denotes the polar cap, and N 
denotes the notch line. 


The main caustic, C : This feature is the result of bunching of emission directions occurring 
due to aberration and time-of-flight effects. 

The secondary (dynamical) caustic, D\ This feature is brought about purely through the 
effects of aberration and time-of-flight, and at higher a merges with the above caustics. This 
feature is accompanied by a spreading of emission directions in the rest of the rotation ((</>,£) 
between (0.3,50°) and (1,0 /dg)\ a trough ). 

The polar cap, P : The polar cap is delimited by the footpoints of the last magnetic field lines 
to close within the light cylinder. Due to the screening of the electric field above the polar cap by 
the particles produced in the pair cascade, no radiation is expected above the polar cap. In LCs 
cutting this feature a very sharp fall in emission intensity can be observed at the appropriate 
phases. 

The overlap region, O: This region is where the distortion of the magnetic field under the 
rotation of the pulsar causes emission directions to overlap. It is associated with the notch in 
the polar cap rim. 

The notch line, N : This line on the phaseplot is a result of the bunching of emission directions 
due to the bunching of the magnetic field lines of the pulsar. 

4. Revisiting the Atlas 

4-1. A single phaseplot 

The LC atlas is produced by generating a phaseplot for each step in the range of possible a , 
and for each of those phaseplots cutting it at each step in the range of possible ( (see Sec. [ 2 ]). 
Fig|3] shows the same phaseplot as is used in Fig. [2j for multiple £ cuts. The colours of the dots 
indicate which feature on the phaseplot the LC peaks are associated with. 

We can now very clearly see the structure of the phaseplot reflected in the set of LCs produced 
from it. We see that at £ < 40° the only dominant peak is associated with the secondary caustic 
(blue dots), except for £ = 0°, where the LC is entirely made up of low level emission. At £ = 50° 
we see that (besides coincidentally cutting through the polar cap (black dots) and seeing the 
accompanying dip in intensity at 0 = 0 and 1) the main caustic (red dots) suddenly becomes 
the dominant feature on the LC as we cut it tangentially at its lowest point. Furthermore, 
we also see the notch line (green dots) for the first time. From £ = 60° onward we see that 
predominantly two-peaked profiles are produced as we are now cutting the caustic above its 
lowest point and our line-of-sight is intersecting with it twice. At £ = 80° and 90° we see pure 
two-peaked profiles. Notice the movement of the peak associated with the notch line moving 
toward earlier phases as we increase £ until its effect becomes lost under the intense emission 
coming from the caustic. Also, notice that the separation between the two dominant peaks (A) 





Figure 3. The set of LCs produced by cutting the same phaseplot (as in Fig. [2]) at successively 
higher £, starting at ( = 0°, and going as high as ( = 90°. Note that the phaseplot is symmetric 
under a 0.5 shift in phase accompanied by a reflection about the equator due to the symmetry 
of the pulsar’s magnetic field (between hemispheres). This means that cuts above ( — 90° are 
redundant. 


steadily increases to almost 0.5 in </> as we progressively cut higher above the lowest point of the 
caustic. 

4.2. The entire set of phaseplots 

Looking at Fig. [4] we see that the intuition we have built up enables us to very easily pin down 
some major attributes of the structure of the atlas which would be produced under a certain 
set of model parameters, in this case a pulsar with a period of 160 ms, a TPC gap width of 5%, 
and a radial gap extent of 1.2i?LC- 



Figure 4. Phaseplots for a between 0° and 90° with a resolution of 30° in both a and ( for the 
TPC model. 


We start with a flat and relatively featureless phaseplot at a = 0° with the structure of the 
magnetic field being the only factor determining the structure of the phaseplot. At a = 30° 
we see the secondary caustics and caustics starting to form, as well as the polar cap. We 
expect to see very few two-peaked profiles, and predominantly profiles with a single early peak 
associated with the secondary caustic. The caustics are much more prominently visible at 
a = 60°, and we can see the notch line and associated overlap region making their appearance. 
The caustics are more extremely curved, and the secondary caustic and caustics have started to 
merge. Here we expect to see a much larger number of two-peaked profiles, as well as some very 
complicated profiles between a = 50° and 80°, having multiple additional small peaks due to the 




























































cut intersecting many of the minor features on the phaseplot. At a — 90° the caustics dominate 
the atlas, yielding a very large range of two-peaked profiles with very large seperations between 
the peaks (A ^ 0.5 in phase). Here the high altitude emissions produce a significant wing 
(around (</>, £) = (0.5,90°)) leading the main caustic at higher £, producing LCs with broader 
peaks. 

5. Discussion and Conclusions 

In this paper we have illustrated how an understanding of the structure of the phaseplots 
generated by geometric pulsar models can translate into an understanding of the structure of the 
resulting LC atlas. Moreover, this understanding also allows us to more easily and effectively 
understand how changes in the model parameters can effect the LC atlas by understanding how 
they effect the phaseplots. This approach also allows us to effectively interpret obtained LC fits 
upon inspection in terms of the features contributing to the best-fit LC, and thus make more 
successful use of the obtained LC fit in, e.g., predicting at which altitudes the relevant radiation 
originates. 
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